rho_phi = dyn_in.M_.params(1);
rho_b = dyn_in.M_.params(2);
rho_th = dyn_in.M_.params(3);
rho_g = dyn_in.M_.params(4);
rho_mu = dyn_in.M_.params(5);
sig_phi = dyn_in.M_.params(6);
sig_b = dyn_in.M_.params(7);
sig_t = dyn_in.M_.params(8);
sig_i = dyn_in.M_.params(9);
sig_g = dyn_in.M_.params(10);
sig_mu = dyn_in.M_.params(11);
slop_p = dyn_in.M_.params(12);
psi_i = dyn_in.M_.params(13);
psi_p = dyn_in.M_.params(14);
psi_g = dyn_in.M_.params(15);
psi_y = dyn_in.M_.params(16);
z = dyn_in.M_.params(17);
phi_p = dyn_in.M_.params(18);
sig = dyn_in.M_.params(19);
psi = dyn_in.M_.params(20);
bet = dyn_in.M_.params(21);
alph = dyn_in.M_.params(22);
del = dyn_in.M_.params(23);
theta = dyn_in.M_.params(24);
PISS = dyn_in.M_.params(25);
MC = dyn_in.M_.params(26);
phi_k = dyn_in.M_.params(27);
phi = dyn_in.M_.params(28);
v = dyn_in.M_.params(29);
b = dyn_in.M_.params(30);
a = dyn_in.M_.params(31);
tau = dyn_in.M_.params(32);
lgv = dyn_in.M_.params(33);
iSS   = z * PISS * lgv/bet ;
RSS   = z*lgv/bet;
MC    = (theta-1)/theta ;
NoK   = ((RSS-1+del)*(1/MC * (1/alph) * (a*b*z)^(alph-1)))^(1/(1-alph)) ;
WSS   = MC * (1-alph)*a^(1-alph)*(z*b)^(-alph)*(NoK)^(-alph) ;
YoK   = (RSS-1+del)*(1/MC*(1/alph)*(1/z)) ;
CoK   = (1-0.2)*YoK-1+(1-del)/z ;
KSS   = (WSS*(1-tau)*phi/(v*b^psi)*(NoK)^(-psi)*(CoK)^(-sig))^(1/(psi+sig)) ;
NSS   = KSS*NoK ;
CSS   = KSS*CoK ;
YSS   = KSS*YoK ;
lamSS = phi*CSS^(-psi)/lgv ;
GSS   = log(0.2*YSS) ;
inveSS = KSS*(1 - (1-del)/z) ;
NSS   = KSS*NoK ;
phi_p = (theta-1) / slop_p ;
xsteady(1) = CSS ;
xsteady(2) = YSS ;
xsteady(3) = GSS ;
xsteady(4) = KSS ;
xsteady(5) = NSS ;
xsteady(6) = WSS ;
xsteady(7) = RSS ;
xsteady(8) = lamSS ;
xsteady(9) = 0 ;
xsteady(10) = 0 ;
xsteady(11) = PISS ;
xsteady(12) = iSS ;
xsteady(13) = MC ;
xsteady(14) = 0 ;
xsteady(15) = lamSS ;
xsteady(16) = RSS ;
xsteady(17) = CSS ;
xsteady(18) = NSS ;
xsteady(19) = WSS ;
xsteady(20) = YSS ;
xsteady(21) = KSS ;
xsteady(22) = log(iSS) ;
xsteady(23) = log(PISS) ;
xsteady(24) = log(z) ;
xsteady(25) = log(z) ;
xsteady(26) = log(z) ;
xsteady(27) = log(z) ;
xsteady(28) = inveSS ;
xsteady(29) = 0 ;
xsteady(30) = 0 ;
xsteady(31) = b*NSS ;
xsteady(32) = log(iSS)-log(PISS) ;
 out(1)=rho_phi;
 out(2)=rho_b;
 out(3)=rho_th;
 out(4)=rho_g;
 out(5)=rho_mu;
 out(6)=sig_phi;
 out(7)=sig_b;
 out(8)=sig_t;
 out(9)=sig_i;
 out(10)=sig_g;
 out(11)=sig_mu;
 out(12)=slop_p;
 out(13)=psi_i;
 out(14)=psi_p;
 out(15)=psi_g;
 out(16)=psi_y;
 out(17)=z;
 out(18)=phi_p;
 out(19)=sig;
 out(20)=psi;
 out(21)=bet;
 out(22)=alph;
 out(23)=del;
 out(24)=theta;
 out(25)=PISS;
 out(26)=MC;
 out(27)=phi_k;
 out(28)=phi;
 out(29)=v;
 out(30)=b;
 out(31)=a;
 out(32)=tau;
 out(33)=lgv;
